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ABSTRACT 

The formation of supermassive black holes is still an outstanding question. In the quasi-star 
scenario, black hole seeds experience an initial super-Eddington growth, that in less than 
a million years may leave a 10 "^ — 10 ® Mq black hole at the centre of a protogalaxy at 
z ~ 20 — 10. Super-Eddington accretion, however, may be accompanied by vigorous mass 
loss that can limit the amount of mass that reaches the black hole. In this paper, we criti¬ 
cally assess the impact of radiative driven winds, launched from the surface of the massive 
envelopes from which the black hole accretes. Solving the full wind equations coupled with 
the hydrostatic structure of the envelope, we find mass outflows with rates between a few 
tens and 10 ^ Mq yr“^, mainl y powered by advection lumi nosity within the outflow. We 
therefore confirm the claim by lOotan, Rossi & Shayjyl (1201 Ih that mass losses can seyerely 
affect the black hole seed early growth within a quasi-star. In particular, seeds with mass 
> 10“* Mq can only form within mass reseryoirs > 10^ Mq, unless they are refilled at huge 
rates (> 100 Mq yr“^). This may imply that only yery massiye halos (> 10® Mq) at those 
redshifts can harbour massiye seeds. Contrary to preyious claims, these winds are expected to 
be relatiyely bright (10^^ — 10^^ erg s“^), blue (Teff 8000 K) objects, that while eluding 
the Hubble Space Telescope, could be obseryed by the James Webb Space Telescope. 
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1 INTRODUCTION 


Supermassive black hole formation is an outstanding question 
in astrophysics. The compelling evidence that links supermas- 


Masorrian et al 19981: 

Ferrarese & Merritt 

2 OOOI: iTremaine et al 

2002; 

Marconi & Hunt 

200.3; Marine & Rix 

2 OO 4 I; Gultekin et al 

2QW; 

McConnell & Ma 

l2013h. stronelv suggests that the answer 


must be sought in the broader context of galaxy assembly. 

A few different formation scenarios have been proposed. Su¬ 
permassive black hole seeds might have a classic al stellar origin 
as the left over of the first generation of stars (e.g. Madau&Reed 
I 2 OOII : IVolonteri. Haardt & Madaull2003l : iTanaka &^aim^ 2009tl 


or they might represent the outcome of the core collapse of 
primordial nuclear star clusters (e.g. (puinlan & Shapirol 


Devecchi_&Wloir^ 120091: Davies, Miller & Bellovard 
2 OI 2 I: ILudI e 


Devecchi et al. 20121 : ILudI et al 2014h . Both proposals 


199C 


2011 


how¬ 
ever, face difficulties in explaining the few observational con¬ 
straints available, namely the early occurrence of massive (> 
10® — 10^® MfVi) accreting black ho les in « ~ 6 — 7 quasars 
i Fan et al1l2006l : lMortlock et al.ll201ll : but see also iTreister et alj 
l201.lh . The main reason is that both scenarios tend to predict 
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small seeds (up to ~ 1000 Mq at most) which are unlikely to 
grow fast enough to power those high-z quasa rs, unless sustained 


super-Eddington accretion is advocated (e.g. 

.Johnson & Bromm 

2007 

; Pelunessv. Di Matteo & Ciardil 20071; 

Milosavlievic et al 

2009 

; Alexander & Nataraianll2014l:lMadau. Haardt & Dottill20l4 

Volonteri. Silk & Dubus||201.'5h. 


Although these observational constraints do not necessar¬ 
ily rule out those mechanisms on a physical base, they may 
more easily be explained by the so called “direct collapse” 


scenario te.g. Bromm & Loebll2003; Begelman. Volonteri & ReesI 

20061; Lodato & NataraianI 

200fil; iDi 

kstra et alJ I 2 OO 8 I; 

Begelman & Shlosman 20091; 

Latif et all 

2 OI. 3 I; Maver et alJ 

201.511. This latter envisages 

a large mass of pristine gas 


10® — 10^ Mq), promptly accumulated at the centre of 
a galaxy-size halo on (sub)parsec scales. A large fraction of it 
would rapidly (< 10® yr) form a massive seed (lO'^ — 10® Mq), 
directly at the centre of a galaxy at z ~ 15. Attractive as it is, 
this process is far from being proved and at least two major steps 
require further investigations. 

Although in principle there is plenty of gas available at 
high redshift and cold flows have been shown to be effec¬ 
tive in bringing tha t gas down to the centre of (massive) ha¬ 
los dpi Matteo et alJ l2012h . the conditions to avoid substantial 
fragmentations and to overcome the centrifugal barrier are not 
fully understood yet. Several possibilities have been discussed. 
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such as the dissociation of H 2 molecules hy Lyman-Wemer ion¬ 
ising radiation coming from nearby, star-f orming galaxies in or¬ 
der to avoid cooling and fragme ntation jperrara & LoebI l2013l : 
iDiikstra. Fenara & Mesingeill2014h . Another possibility is the on¬ 
set of supersonic turbulence and the removal of angular mo¬ 
mentum due to non-axisymmetric perturb ations and gravitational 
torques during the collapse of the halo feegelman & ShlosmanI 


l2009l : lchoi. Shlosman & Begelmanl2013l . l201^ . or at the centre of 

major m erger remnants between rare and massive galaxies at high 
redshift jMaver et aDl20ld.l20Bl . 

The second issue is how to actually form a black hole and what 
is its initial mass. The answer may vary according to the physical 
properties of the assembled mass. When more than ~ 10® Mq 
can be rapidly piled up, the resulting structure likely becomes dy¬ 
namically unstable (even if rotating) and relativistic radial i nstabil¬ 
ity can lead to implosion and direct black hole formation jpowleil 
ll966l : lBaumgarte & Shapircll999l : IShibata & ShaDiroll2002l) . How¬ 

ever, forming such a structure requires r ather extreme conditions 
(e.g. a major merger: [Mayer et ^l2015h . With relatively milder 
accretion rates, nuclear burning can start at th e centre of a con - 
vectively stable object, i.e a supermassive star llBegelmaiill201C)h . 
After a million years, the core that could not convectively ac¬ 
quire fresh hydrogen collapses to form a stellar size (~ 100 Mq) 
black hole. Highly optically thick gas keeps however falling onto 
the newly born black hole, with enough angular momentum to be 
able to generate accretion power. This energy feedback inflates 
the innermost part of this inflow, creating a quasi-star, a mas¬ 
sive, slowly rotating envelope, sustained against its own gravity 


by th e black hole accre t ion power (|Begelman,_Rr3ssj_&_Armitage 


200S; BegelmaiJ 2010l: Volonteri & BegelmanI 201fll : 


__ __ _ __ iBall et al, 

TOlltlDotamRossT^Shavi^lTOnlk At ouasi-star centres, the em¬ 


bryo black holes may accrete at a super-Eddington rate, as en¬ 
ergy is transported outward by convection (not by radiative dif¬ 
fusion) through the envelope. The initial expectation was that 
seeds of ~ 10** — 10® Mfy, ma y easily grow in < 1 Myr 
jBegelman. Rossi & Armitaral2008h . After this time, the envelope 
would be definitively dispersed and accretion would proceed at an 
Eddington limited fashion directly from the protogalactic d i sc. 

This result was questioned bv IPotan. Rossi & ShavivI (1201 ih . 
that tried to quantify the impact of radiative driven winds from the 
surface of these weakly bound envelopes. They found that in a large 
part of the black hole mass-envelope mass parameter space, winds 
can be so powerful that the envelope evaporates before the black 
hole is able to double its mass. This effect would greatly limit the 
number of protogalaxies in whose centre the conditions are prone 
to massive > 10"^ Mq seed formation. Moreover, because in their 
model most of the radiation energy in diffusive luminosity is con¬ 
verted into kinetic energy of the wind, quasi-stars would be very 
dim objects, virtually undetectable. In that paper, however, the ra¬ 
diative driven wind was not modelled solving the full equations of 
motion and in particular, the advection energy term was neglected. 

We therefore set out to critically reconsider continuum driven 
winds from the surface of radiation-dominated objects. In fact, we 
find that the advection term has an important dynamical role as the 
main driver of the wind. We then explicitly consider quasi-stars 
and calculate mass loss rates and photospheric properties and we 
assess their detectability with Hubble and the lames Webb Space 
Telescopes (IWST). Our result is that although winds are still a 
major limit for black hole growth, photospheric luminosities in the 
wind ensure possible bright targets (10''"‘ —10'*® erg s“*) for IWST. 

This paper is organised as follows: in Section [2 we describe 
our model of radiation-dominated wind and we discuss the main 


properties; in Section[3]we couple this wind prescription to the hy¬ 
drostatic envelope predicted for quasi-stars, finding equilibrium so¬ 
lution and discussing their evolution. Section|4]is devoted to simple 
predictions regarding the observability of quasi-stars by current and 
future space-based telescopes. We discuss and summarise our main 
findings in Section [2 addressing the main limitations of our work 
and future steps. 


2 THE WIND MODEL 

2.1 Equations and general properties 


We consider a stationary, spherically-symmetric, radiation- 
dominated wind launched by a non-rotating, stellar-like object of 
mass M* from a spherical surface of radius i?*, that represents 
the base of the wincQ By radiation-dominated, we mean that the 
contribution of the gas pressure pgas is assumed to be negligible 
compared to the radiation pressure PraA, i-e- Pgas/Prad 1. This 
assumption allows us to neglect the presence of pgas in the follow¬ 
ing calculations (and we simply write p = Prad), but it requires at 
the same time that the wind is launched from a radiation-dominated 
object (as we consistently show for quasi-stars in Section|2l- Just 
outside Rt, the gas is assumed to be initially optically thick and in¬ 
teracts with radiation through a constant opacity k. In this Section 
we implicitly assume the Thompson scattering opacity (though the 
specific value is in fact irrelevant for our results since we develop 
all our calculations in a dimensionless form), while in Section|2we 
will adopt a temperature-dependent opacity law (see Section lrTT l. 
We stress, however, that the (dynamical) results presented in this 
Section are largely independent of the assumed opacity law, as also 
discussed in Section [33] We are interested in primordial composi¬ 
tion objects, where line-driven interaction is negligible. The equa¬ 
tions that describe this system are similar to those used by several 
previous works about stellar winds (and spherical ac cretion) in both 


the optically thin and optically thick regime (e.e. Zvtkow 

I 1972 I 

lBegelmanlll978lll979l:lKatolll983l:IOuinn & Paczvnskilll985 

): 


M = 4nFp V, (1) 

1 GM^ kL 

2 dr ^2 


mA (t. _ AA + p±!l\ = _ A 

dr \ 2 r P ) dr 


(3) 


These equations determine the structure of the gas density p, the 
radial gas velocity «, the luminosity carried by photons L, the (ra¬ 
diation) pressure p and the (radiation) internal energy density U as 
a function of the spherical radius r within the gravitation poten¬ 
tial d? = —GMi^jr induced by M* outside Ji*. The steady-state 
wind is characterised by the constant outflow rate M. Equations]!] 
[2and|2describe the conservation of mass, momentum, and energy, 
respectively. 

Such a system of equations is not closed and several ap- 
proach es can be used to cl ose it to different degrees of app roxi¬ 
mation. IshavivI (l2001bh and lOwocki. Gavlev & ShavivI (l2004l start 
from similar equations, except that they initially include the con¬ 
tribution of the gas thermal pressure to the momentum and energy 
conservation. Then, they simplify the system focusing on the super¬ 
sonic branch, thus subsequently neglecting the gas pressure terms 


^ In the following, we will always use the subscript * to indicate quantities 
evaluated at fj* 
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Radiation-dominated winds and quasi-stars 3 


in the momentum equation (which brings it back to our same equa- 
tion|3 and the advection term (p + U)/p (with the corresponding 
one due to gas pressure) in the energy equation. This approxima¬ 
tion leads to the great advantage that fully analytic solutions can 
be derived. However, the limitation is that the lack of the advec- 
tive term makes the behaviour of the wind insensitive to the local 
optical thickness. 

Instead, we fo ll ow an approach similar to that used by e.g. 
lOuinn & Paczvns^ ill 98.511 and we explicitly include additional 
prescriptions to properly describe the behaviour of the wind in the 
extrema of very optically thin and optically thick regime, i.e. when 
the optical depth: 

^ + oo 

r(r) = / AC p(x) dx, (4) 

J r 

is either r <C 1 or r 1, respectively. When the wind is opti¬ 
cally thick, radiation and matter can reach local thermodynamical 
equilibrium at the same temperature T (which relates to the energy 
density U = 3p — aT'^, where a is the radiation constant) and the 
gradient of the radiation energy density is: 


dr 


T:g>i 


3KpL 

dTrr^c 


(5) 


On the other hand, local thermodynamical equilibrium may not be 
reached in the optically thin limit and a unique temperature may not 
be a physically-motivated quantity. In this case, the photons carry¬ 
ing L travel with roughly radial orbits and i nteract very lit tle with 
matter, keeping L almost constant (see e.g. IZvtkowl[l972h . Then, 
the radiation energy density decreases mostly because of geometri¬ 
cal dilution in a progressively larger volume: 


dr 


( 6 ) 


T<1 


27rr3c 

We follow lOuinn & Paczvnskil (Il98.5h in defining the total gradient 
of U as the sum of the two limiting cases: 

^ L 
dr 


/(t 


2Tir^c 

where we define the function: 
St” 

= — + f = iipr. 


(7) 


(8) 


The “effective” opacity r leads the gradient of U to the right opti¬ 
cally thin and optically thick limits when r <C 1 and f ^ 1, re¬ 
spectively. However, f is just an approximation of the actual opac¬ 
ity r; the two are related by a constant factor when r is a power law 
and such a factor is close to 1 when r oc r~^. Although we do not 
know a priory the relationship between r and f, we demonstrate in 
the following that r oc r~^ roughly holds and therefore r ~ r (see 
Sectionl^. 

Finally, we need to re l ate p a nd U to close the system of equa¬ 
tions. lQuinn_&^aczyn^ lll985h implicitly assume that p = Ul3 
everywhere in the flow (see their equations 11b and 12). This is 
correct in the optically thick regime, but is not valid when the 
gas is optically thin. Indeed, U — p when the gas is optically 
thin; this different relation between p and U is also responsible 
for the inexact relation between the luminosity observed by an ob¬ 
server at in finity and by an obs e rver c omoving with the flow, as 
reported by ^uinn & Paczvnsl^ d 198.51 : see also Section 12.21 and 
ICassinelli & Castoifig?^ In order to have a smooth transition be¬ 
tween the two regimes, similar to the case of the gradient of U 
(see equation|7j, we propose the following functional form for the 
opacity-dependent ratio p/U'. 



Figure 1. The tick continuous and dashed lines show the behaviour of / 
and g as a function of t, respectively. For reference, the vertical dotted and 
dashed lines mark the points f = 1 and f = 2/3, respectively. 

Figure [T] shows the behaviour of both /(r) and (/(r). The latter 
goes correctly from 1 (when the wind is optically thin) to 1/3 (when 
the wind is optically thick). However, the accuracy of both pre¬ 
scriptions is questionable around f ~ 1, because the actual form 
is largely arbitrary. We therefore compute and compare wind mod¬ 
els, choosing different functional form for /(r) and g{f) and con¬ 
cluded that our results are not affected as long as the the limits are 
correct and the transition occurs rapidly (over several r) around 
f ~ 1. 

The equation of the conservation of energy can be directly in¬ 
tegrated, becoming an algebraic equation for L: 

L{r) ^ E - M ^ + {1 + g{f))^y (10) 

where we make use of g{f) explicitly and we introduce the total 
conserved luminosity E, which represents the constant of integra¬ 
tion. The system of equations that we finally solve (often dubbed 
as “wind equations” in the following) is composed of equations [T] 
[2|7]and[T0l coupled with the definitions of f, /(r) and gij). The 
system has two ordinary differential equations and two algebraic 
equations for the dependent variables p (or r), v, L and [/ as a 
function of r. 

At this point, it is convenient to introduce new dimensionless 
variables. We define the new velocity variable w = w^/Vesc^ where 
Wesc = 2GMt /Rt is the escape velocity from the base of the wind, 
the new radiation energy density variable u = UkR^ / (GMi,), and 
the new luminosity variable T = L/L^^^. F is the Eddington ratio 
and the Eddington luminosity Z/Edd is defined as: 

LEdd = = 1.26 X 10^® m* erg s~^, (11) 

K 

where k is the opacity in units of the electron scattering opacity 
Kes ~ 0.35 cm^ g“^ (assuming primordial abundances) and m* 
is the stellar mass M* in units of solar masses. The independent 
variable r can also be transformed into x = 1 — i?*/r, such that 
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4 Fiacconi & Rossi 


the interval r G [R-t, +oo) is mapped into x G [0,1). We can first 
express equation]^ with the new variables as: 

= ( 12 ) 

where here and in the following ' — d/dx. The gradient of the 
radiation energy density u becomes: 

«'=-2 (1 - x) r/(f), (13) 

where /(f) is defined in equation[8]and: 

(14) 


procedure outline d bv lOuinn & Paczvnskil ( Il985l : see also IZvtkowl 
ll972l : [K^ll983h . First of all, we characterise our star by choosing 
a value for /3. Then, we pick a value for a. and £. These two con¬ 
stants of integration are not enough to fully characterise the wind. 
We need a boundary condition, specifically the asymptotic gas ve¬ 
locity at infinity Wao- With that, we can specify the initial condi¬ 
tions at a large radii x^a = 1 — <5, much larger than i?* (i.e. when 
(5 —>■ 0) and start our integration of the wind from outside inward. 
Practically, we calculate the radiative luminosity as seen by an ob¬ 
server at infinity as: 


CjdO — ^ CXIVo 


(18) 


is the definition of f using our dimensionless variables. We intro¬ 
duce the two factors a and /3; a is a d imensionless expressio n for 
M in terms of the trapping radius Ru llBegelmanll 1978111979h : 


Rti kM 
Ri, TtfcT?* 


(15) 


The trapping radius is where the diffusion time scale for photons is 
equal to the dynamical time of the outflowing wind, i?tr(f/c) « 
Rti/v, which implies that f « c/v at Rti. We will see in the fol¬ 
lowing that within this radius, since the radiation is trapped, the 
luminosity transported by diffusion becomes subdominant with re¬ 
spect to the energy advected within the flow. The parameter /3 is a 
dimensionless factor depending on the properties of the star only: 

/3 = —^ « 486 (16) 

'^^esc 


where r* is the stellar radius in units of solar radii. P appears natu¬ 
rally from the normalisation of the optical depth and measures the 
deepness of the gravitational potential well of the parent star. It de¬ 
pends on the stellar properties 77* and M* only and therefore, M* 
and P (or Uesc) are enough to rescale the equations in physical units. 
Finally, the algebraic equation for F can be obtained from equation 
[T^and reads: 




-\- X — 1 -\- 


1 +g(f) 

1 — a: f / ’ 


(17) 


where S = E/L-Edd and gij) is defined in equation[9l The depen¬ 
dent variables w, u and F are proportional to the kinetic energy of 
the gas, to the energy density of the radiation and to the luminosity 
carried by photons, respectively. Therefore, equations ! 121113l andl 17 1 
compose the system that describes the energy exchanges between 
the different components of the system. 


2.2 Numerical integration of the wind equations 

We integrate numerically the wind equations (equations! El [13 


and M. using the CVODE mo dule of the SUNDIAL 


ions [2] 
lAuj^ p; 


age dCohen, Hindmarsh & Dub^ 1 19961 : iHindmarsh et al 


pack - 


20o3) . 


CVODE is a C solver for stiff and non-stiff ordinary differential 
equation systems in explicit form. We adopt a fifth-order back¬ 
ward differentiation formula in fixed-leading coefficient form with 
a modified Newton iteration to solve non-linear systems. CVODE 
provides also a module to find the roots of nonlinear equations 
which is well suited to determine the position of the photosphere 
and the local properties of the wind while contemporary solving 
the wind equations. 

In order to find solutions of the wind equations, we follow the 


At Xoo, the luminosity Foo comoving with the flow can be consid¬ 
ered constant and given by: 


Foe = 


£o 


1 -F 2wli^/P 


£ — awoo 

l + 2w]i'^/p' 


(19) 


We can then integrate the radiation energy density in the optically 
thin limit: 


w'~-2(1 - a;)Foo ^ M = Foo(l-2:)^ (20) 


where we use the boundary condition m( 1) = 0, and the wind ve¬ 
locity: 

Ml'= Foo — 1 => u; = Woo — (Foe — 1)(1 — a;). (21) 


We can also write the explicit behaviour of r: 


aP{l — x) 

\/Woa - (Foe - 1)(1 - X) 


( 22 ) 


Note that equation [Incomes out naturally by evaluating equation 
I17lat a; —> 1, where 1 + g{j) —>■ 2 and u and r are described by the 
expressions above (equations 120! and [2^ . Next, we choose a value 
for 5, typically <5 ~ 10“®, we check that indeed t(xoo) <C 1 and 
we use the formulas above to provide the initial conditions w(xoo) 
and u(xoo) for the wind equations. 

We then integrate the equations inward up to the surface x = 0 
(or up to the point where a solution exists). We define the photo¬ 
sphere as the place where the equality Lphot = 4:TrRpEotttTpEot 
is satisfied, where a = ca/4 is the Stephan-Boltzmann constant, 
whereas the temperature T is defined from the energy density U as 
T = ([//a)^'^'*, regardless of the local optical depth. Such a tem¬ 
perature is a proxy for the local temperature and it recovers its full 
physical meaning only when f > 1. The photosphere identified in 
this way usually lays at f ~ 2 — 3. 

Every solution of the wind equations is specified by the pa¬ 
rameters a, £ and Woo, once the underlying star is set by p. Among 
those parameters, a is directly related to the outflow rate M and is 
necessary to solve the wind equations, i.e. such a model does not 
allow to infer theoretically the value of M a priori. However, an ac¬ 
ceptable solution has to satisfy additional self-consistency require¬ 
ments, which in turns impose constraints of the parameter space 
and ultimately on the value of M. Those self-consistency require¬ 
ments are imposed by the assumption that the wind originates from 
a star-like object. In particular: (i) the solution has to extend in¬ 
ward to at least x = 0 ; (ii) the wind has to be optically thick close 
to the surface of the star, i.e. the photosphere has to be above the 
base of the wind, namely r* > 1 and 77phot > 77*; and (iii) the 
wind has to connect to an hydrostatic solution, i.e. it should be ini¬ 
tially subsonic (i.e. A4* < 1 ) and with a moderate velocit 50 (i.e. 


^ SUNDIALS is publicly available at https : //computation. llnl. 

gov/casc/sundials/main. html. ® As consequence of the assumption of steady-state, we note that we cannot 
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Figure 2. Relation between the velocity at infinity Woo and the Mach num¬ 
ber at the surface A^* assuming /3 = 100, a = 1 and £ = 5. The dotted 
thin line marks Ad* = 0.5 for reference. The range of possible Wao is very 
nan'ow. The inset shows the distribution of the Mach numbers Ad* obtained 
by 500 realisations. A clear peak around 0.8 — 0.9 is present. 


w* <C 1). A wind solution is then accepted only when it satisfies 
all the conditions listed above, and it is discarded otherwise. 


2.2.1 The wind velocity at infinity 

To try and simplify further our procedure, we first asses the sen¬ 
sitivity of our solution to our choice of Woo. The arbitrariness of 
Woo is simply a consequence of our neglecting gas pressure, in 
the equations describing a radiation dominated wind. When ther¬ 
mal gas pressure is explicitly accounted for, the sonic point of a 
solution (i.e. where Ad = 1) is also a critical point (i.e. a diver¬ 
gent point for w'). The requirement on the position of the critical 
point to cure the local divergency translates naturally into a condi¬ 
tion that fixes the value of Woo • As a consequence, solutions of the 
wind equations with gas p ressure only dependent on a and S (e.g. 
lOuinn & Paczvnskilll985h . In our solutions, instead, the sonic point 
is not a critical point, and Woo is not univocally determined. How¬ 
ever, the fact that with gas pressure terms there is only a single value 
for Woo and that solutions should be continuous as Pgas/Prad —>■ 0 
suggests that the range of possible Woo may be nan'ow. Therefore, 
we investigate this possibility. 

We setup a grid of five representative values for /3 € 

{10, 50,100, 500, 1000}, and five representative values for £ G 
(1.5, 2.5, 5, 7.5, 10}. For each pair (/d, £), we divide the interval 
log a G [—3,1] uniformly, and for each value of a we run 10^ in¬ 
tegrations of the wind equations choosing a random value for Woo 
distributed uniformly in the logarithmic interval [—3, log(f/a)]. 
We keep only the acceptable solutions according to Section lZ^ 

Our results confirm that the range of Woo that leads to self- 
consistent solutions is narrow, usually < 0.1 dex, and centred 
around Woo ~ 1- The values of Woo also correlate with Al* in 
the interval 0.1 < Af* < 1. Such a correlation is shown for 500 


allow for m* = 0 because the density would otherwise diverge, as implied 
by the conservation of mass in equation[T] 


■to 


■to 



No winds 

'lQ-3 ' ' .. 

a 


-- 0.31 

1 - 0.6 
-0.9 


Figure 3. Monte Carlo sampling of the parameter space a—f. Upper-panel: 
blue circles, red squares, green triangles and magenta stars show the results 
for 0 = 1000, 100, 10, and 3, respectively. For visualisation purposes, we 
show 10% of the 10^ calculations performed for each value of 0. Lower 
panel: the same as above showing all the results for the case 0 = 100 
where we colour-coded the points according to the ratio between the ad- 
vected luminosity at the surface Ladv * and the radiative luminosity L* 
at R^,. Both panels show the limits to the parameter space imposed by the 
self-consistency conditions as described in the text. 


realisations in Figure|2]for an example configuration with 0 = 100, 
a = 1 and £ = 5 and exhibits typical features common to all the 
other combinations of parameters. In particular, most of the interval 
of allowed Woo corresponds to values of the Mach number larger 
than ~ 0.5 — 0.6, as shown by the distribution of AI* represented 
in the inset of Figure|2 peaking around Al* ~ 0.8. Motivated by 
that, we can use this occurrence as an approximate additional con¬ 
straint to remove the freedom of choosing Woo by choosing a value 
for A4* to be matched at 7?*. In the following, we focus only on 
solutions with Al* = 0.8 ± 0.05. 


2.3 Results 

We are now in the position to explore the properties of the solu¬ 
tions within the parameter space (a, £) as a function of 0. Figure 
[^highlights the boundaries on the parameter space imposed by the 
self-consistency conditions. Solutions are limited from below by 
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6 Fiacconi & Rossi 


the requirement that £ ^ 1. We assume this condition as neces¬ 
sary in order to launch the wind. In fact, the actual condition to 
have an accelerating wind is F* > 1 or in other words, that the star 
should shine above the Eddington limit (equation WH . Nonethe¬ 
less, we conveniently chose the limit £ ^ 1 because (i) it implies 
r* > 1 and (ii) the hydrostatic solution (the star) below the wind 
naturally provides « L* as a boundary condition (we discuss 
this point with more details in Section]^. The limit at small a and 
F < 10 results from imposing f* > 1 and it depends on fi since the 
normalisation of the optical depth is oc a/3, as shown by equation 
[Tdl and in the upper panel of Figure |3 Physically, this is because 
matter needs to be faster to escape from a more compact star, and 
from mass conservation (p oc M/v) it follows that a higher mass 
loss rate is required to maintain the same optical depth f* > 1. Fi¬ 
nally, at fixed a, the upper value of £ is constrained by the matching 
with a hydrostatic solution below the wind, that requires < 1. 
Incidentally, we note here that our set of prescriptions do not set an 
upper limit on E. This will be provided by the physical characteris¬ 
tics of the stellar object powering the wind, more explicitly by how 
much super Eddington its emission can be. We will explicitly show 
this in Section[3for quasi-stars. 

The maximum mass loss rate for a given £ (grey area on the 
right) is instead a physical limit. This is obtained when the dom¬ 
inant energy source for the kinetic luminosity of the wind is the 
enthalpy of the gas and £ ~ Woo a. In practise, this is the be¬ 
haviour of an adiabatic wind. We show this by considering a fully 
adiabatic solution. This latter has no radiative luminosity in its gov¬ 
erning equations (see Appendix |3, and the presence of a critical 
point allows us to relate the velocity at infinity tUoo.adiab to the 
condition at the base of the wind. We can therefore derive that an 
adiabatic wind with AT* = 0.8 will have Woo.adiab = 3wc ~ 0.6, 
where Wc is the velocity at the critical point Wc = Sc, related to 
AI* by equation I A4I and the relation plotted in Figure IaTI For a 
given £, we therefore have; 

amax«F/0.6. (23) 

This relation is plotted as a dashed line in Figure [3] and clearly 
marks the rightmost limit of our solutions and the beginning of the 
“forbidden” region dubbed “adiabatic limit”. 

That our solutions tend towards an adiabatic behaviour as M 
increases, can be better appreciated by looking at the lower panel 
of Figure[3 which shows solutions for /3 = 100, colour coded ac¬ 
cording to the ratio Fadv,*/F*, where Ladv = M(p + U) / p is 
the luminosity advected within the bulk outflow. Across a ~ 1 
the regime of the wind changes: for a < 1 energy in radiation 
is mainly transported by diffusion while for a > 1 advection be¬ 
comes more and more dominant as the outflow rate increases to¬ 
wards the adiabatic limit. As mentioned before, this is exactly the 
physical meaning of the trapping radius (see also Figure [4] lower 
rightmost panel) and a > 1 implies that the trapping radius occurs 
within the outflow, Rti > R* (see equation fTSt . Formally, a fully 
adiabatic solution has i?tr = oo (i.e. a = oo) and correspondently 
Fadv,*/F* = 00. Advection is indeed the only transport mecha¬ 
nism in an adiabatic wind. We can go a step further and calculate 
M of a solution relative to the maximum possible mass loss rate. 
The latter corresponds to that of an adiabatic wind (Mad, equation 
lATl launched from the same star, with the same initial conditions 
(i.e. the same Cs,* and p* at f?*) at the same £ (which for an adia¬ 
batic windreadsF = Mad(uJ/2-F SCg^* — G'M*/i?*)/I/Bdd): 


M 

Mad 



(24) 


When the contribution of the radiative luminosity L* at i?* to the 
total energetic budget becomes negligible, M approaches the adia¬ 
batic value. 

The upper panel of Figure [3 shows the comparison of 10^ re¬ 
alisations for four different values of /3. We find that the winds orig¬ 
inating from more compact stars (i.e. those with smaller P) sustain 
higher outflow rates at the same total luminosity £. This is again 
related to the fact that gas needs to be faster to escape from a more 
compact star, and it compensates this increase in velocity by an in¬ 
crease in M to maintain optical thickness. As a consequence, an 
optically-thick wind needs to carry a larger and larger fraction of 
the total luminosity in kinetic and advective form. This explains 
why the range of possible M gets squeezed towards the adiabatic 

limit foi0/3 1. 

Figure |4] shows two example solutions of the wind equations 
with P = 100. The upper row shows a solution with a = 0.4 and 
F = 5, i.e. at the centre of the allowed region in the lower panel 
of Figure[3 while the bottom row shows a solution with a = 7.64 
and F = 5, i.e. very close to the adiabatic limit. The wind velocity 
of the first solution grows steeply within the photosphere and then 
flattens to its asymptotic values; on the contrary, the sound speed Cs 
decreases quickly within the photosphere, matching the gas veloc¬ 
ity at the sonic point very close to i?*, beyond which the wind be¬ 
comes highly supersonic. Note that, regardless of the optical depth, 
we always use the general definition of sound speed: 


(^\ ^ If 

\dp)s 3p’ 


(25) 


where S is the specific entropy and the second equality is based on 
the most general equation of state for a non-isentropic, radiation- 
pressure dominated fluid, namely p{p,S) = K{S)p^^^. In fact, 
this relation is only valid when p = (7/3, i.e. in the optically-thick 
limit, but we need to consider the behaviour of Ca only within the 
photosphere. The solution close to the adiabatic limit shows a simi¬ 
lar behaviour, though the gas velocity grows less steeply than in the 
previous case and mostly within Ru. It also reaches an asymptotic 
velocity lower than the previous case, though the total energy £ is 
the same. 

The profiles of the optical depth, density and radiation energy 
density are similar in both examples. We can fairly accurately de¬ 
scribe them with power laws, at least close to and outside the pho¬ 
tosphere. In particular, the optical depth decreases with radius as 
r~^, while the density as r~^ once the gas velocity remains almost 
constant. This confirms a posteriori that r is a good approximation 
for T. We note also that the wind is highly optically thick close to 
i?*, while the photosphere forms at Tphot ~ 2 — 3. This latter con¬ 
firms that our approximate treatment of the transition between the 
optically thin and thick regimes gives sensible results, close to the 
conventional Tphot ~ 1- 

The largest difference between the two solutions is in the en¬ 
ergy budget, shown in the left column of Figure|4] There, we com¬ 
pare the profiles of the radiative luminosity L, the advected lumi¬ 
nosity Ladv, and the kinetic luminosity Lkin = Mv^ {2, all nor¬ 
malised to Z/Edd for convenience. These luminosities, summed up 
with —GMMi,/r, give the total, constant luminosity E. The en¬ 
ergy budget in the first solution is dominated by the radiative lumi¬ 
nosity from the base of the wind (indeed Rti < 7?*) to infinity. Af- 


Formally, such a limit cannot be reliably modelled by our wind equations, 
because it would imply approaching a regime where general relativistic cor¬ 
rections might become relevant. 
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Figure 4. Example solutions of the wind equations for 0 = 100. The upper row shows the solution for a = 0.4 and £ = 5, while the bottom row shows 
the solution for a = 7.64 and €, very close to the adiabatic limit of Figure[3] Left column: continuous, dashed and dotted lines show the profile of v/v^bc, 
Cs/uesc and M, respectively. Central column: continuous, dashed and dotted lines show the profile of f, u and pjp-k, respectively. Right column: continuous, 
dashed and dotted lines show the profile of ikin/i'Bdd. ^adv/^Edd and respectively. The dark and light grey shaded regions in all panels mark 

the base of the wind and the surface of the photosphere, respectively, while the red, vertical, dash-dotted lines in the bottom row indicate the position of the 
trapping radius. The thin lines in the bottom-left panel and the thin dotted line in the bottom-central panel show the velocity, sound speed, Mach number and 
density, respectively, of the adiabatic wind with the same critical point of the wind solution. The two solutions are very similar within fJtr. 


ter a small decrease, it remains almost constant with radius, while 
ikin quickly rises within the photosphere, yet remaining subdom¬ 
inant. This behaviour indicates that the acceleration of the wind 
is not powered by L but it rather occurs at the expenses of Ladv, 
that drops accordingly within the photosphere. In contrast, Ladv 
in the a = 7.64 case (lower panel) is initially an order of mag¬ 
nitude higher than both L and Lkin and becomes comparable to 
L around J?tr. Instead, L remains always subdominant compared 
to Z/kin- The radiative luminosity’s mild growth within the photo¬ 
sphere is similar to that of an “effective” Eddington ratio that one 
would obtain from equation[T2l Tefi = 1 + w', when w' from the 
adiabatic equations IaT] and lA2l is used. Notably, also in this case, 
Lkin is ultimately powered by L^dv 

Finally, we show explicitly the similarity between the solution 
close to the adiabatic limit and an actual adiabatic solution in the 
bottom row of Figure |4] Specihcally, after calculating the solution 
of our wind equations, we calculate also the adiabatic solution (see 
Appendix lAt which has the same critical radius. Such an adiabatic 
solution crosses naturally the critical point with the same value of 
critical velocity Wc of the full solution, and closely resembles it 
(comparing v, Cs, M, and p) within the trapping radius, outside of 
which radiative diffusion becomes relevant. 

From our results, it is clear that the ultimate source of kinetic 
energy for the gas is advection energy and not the diffusive radi¬ 
ation luminosity. This is not consumed as the gas propagates out¬ 
ward. To increase the mass loss, it is therefore necessary to go to¬ 
wards an adiabatic solution, where initially the advection energy 


domina t es the energy budget. This is different from the results by 
IShavivI ll2001bli . [Owocki,_Ga^le^^^haTO (2004), and from the 
description used in iDotan. Rossi & ShavivI ( 120 llh of a “photon- 
tired” wind. 


3 WINDS FROM QUASI-STARS 


Directly forming massive black hole seeds is possible in princi¬ 
ple when a mass of ~ 10® Mq can be collected in ~ 10® yr, 
requiring inflow rate > 100 Mq yr“^. The reason is that such a 
rapid accumulation of mass has to occur before nuclea r reactions 
domi nate the evolution, setting a lifetime of ~ 10® yr teegelmaiil 
l2010h . In these extreme conditions, a dynamical instability due 
to relativistic effects may develop and even rotation can not pre¬ 
vent t he cloud from collapsing dir e ctly into a black ho le i Fowled 
I1966L iBaumgarte^^Shagi^ Il999l : IShibata & Shapirol l2002l : see 
howevei lFemrart‘aLll201^ On the other hand, milder conditions 
(e.g. inflow rate ~ 0.1 — 1 Mq yr~^) would lead to the forma¬ 
tion of a supermas sive star, possibly ~ 10® — 10^ Mm, stabilised 


by some rotation (lBegelrnM l201o| : iHosokawa. Omukai & Yor^ 
I 2 OI 2 I: iHosokawa et aljpOlSh . After ~ 10® yr, a small embryo 
seed (< 100 Mq) can form at the centre o f a supermassiv e 
star at the end of the hydrogen-burning phase llBegelmanll201C)h . 
Such a seed then needs to go through a phase of vigorous super- 
Eddington accretion to reach 10^ — 10® Mq within a few mil¬ 
lion years from its birth. This super-Eddington accretion can oc¬ 
cur within a quasi-star: a very massive (> 10® — 10® Mq) quasi 
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hydrostatic envelope that surrounds the black hole and feeds it 
at a r ate equal to roughly its own (i.e. quasi-star’s) Eddington 
limit teegelman. Rossi & Armitagal20Q8h . Accretion around the 
Eddington limit involves radiation dominated gas, which is loosely 
bound with a total energy close to zero. For this reason outflows can 
easily form. Flere we investigate whether this vigorous accretion in 
quasi-stars is also accompanied by mass l oss, as expected in other 
super -Eddington systems such as discs l lBlandford & BegelmanI 
l2004h . 


3.1 Equations 

We follow Begelman. Rossi & Armita^ ( l2008h and 
IPotan. Rossi & Shaviv i201lh to describe the hydrostatic en¬ 
velope of a quasi-star and then we match it with our wind model, 
looking for equilibrium solutions. A quasi-star is made of four 
components: (i) the central, accreting black hole, (ii) a convective, 
radiation-pressure dominated envelope, (iii) a porous radiative 
layer, and (iv) a wind. In th e following, we briefly desc ribe 
components i-iii, remanding to IPotan. Rossi & ShavivI (1201 lli for 
additional details. We do not model explicitly the inflow on to the 
central black hole; instead, we treat it as a boundary condition, 
assuming that a black hole o f mass M, is accreting through a 
convection-dominated disc dStone. Pringle & BegelmanI 199S; 


]^umenshc hev & Abramowicj 199^ OuataerT&Gmzin^ 20QC ; 

Agol et all 200 ill within a few Bondi radii rs = GM,/(2Cs^c). 
where Cs,c is the central sound speed < tb. We assume that the 
black hole is radiating at a luminosity L, close to the maximum 
that convection-dominated accretion flows can sustain evaluated at 
5rB, namely: 


Lm — Fconv(5rB) — 4 TT (SfB) pc 


(26) 


where pc is the central density outside a few tb. L, it is injected 
at the centre of the envelope and transported till the base of the 
wind, first by convection and then by diffusion. We used it as an in¬ 
ner boundary condition for the integration of the envelope and we 
also neglect the gas mass within 5rB, assuming that only M, con¬ 
tributes at smaller radii. We checked the effect of varying the po¬ 
sition of the inner boundary. We find that most of the properties in 
the M, — Mi, p lane (see Sectionl3.3b remain unchanged. However , 
as discussed bv iBall et al.l ( 1201 ll) and iBalL Tout & Zvtkowl ( l2012h . 
changing the inner radius modifies (at the same central pressure) 
the ratio M./M* of two consistent solutions. This slightly dis¬ 
places the no-hydrostatic-solution region. Although a few Bondi 
radii are a reasonable estimate for the inner accretion region where 
the black hole gravity is expected to dominate, we caution that such 
a choice remain somewhat arbitrary. 

Outside 5rb, the radiation-pressure dominated, convective en¬ 
velope extends; we assume that it satisfies the hydrostatic equilib¬ 
rium: 

dP GM{r)p 

dr r2 ’ ^ ’ 

where the total pressure P is the sum of the gas pressure Pg = 
pkBT/(pmp) and the radiation pressure Pr = aT'^fS, specified by 
the gas density p and temperature T. ks, a and rrip are the Boltz¬ 
mann constant, the radiation constant and the mass of the proton, 
respectively. We assume the mean molecular weight p = 0.59, ap¬ 
propriate for gas with primordial composition at P > 10^ K, as 
usually true everywhere in the interior of quasi-stars. The enclosed 
mass M{r) is given by: 


M{r) = M,+4TT / p{r')r'^dr'. (28) 

JStb 

The dominant energy transport mechanism within the envelope is 
convection, which induces a temperature gradient very close to adi¬ 
abatic; therefore, we evolve the temperature gradient assuming that 
it is equal to the adiabatic one: 


dlogP _ 7 ad - 1 
d log P 7ad 


where the adiabatic index 7ad depends on the ratio 
according to: 

_ 32 - 24C - 3C^ 

“ 24 - 18C - 3(2 ' 


PJP 

(30) 


The convective envelope extends till the radius Tconv where 
Tconv(fcon) = L,. Outside Tconv, T/conv < L, and diffusion be¬ 
comes more efficient in transporting L, than convection. However, 
L, may be larger than the local Eddington limit associated to the lo¬ 
cal enclosed mass M(r) and to electron scattering opacity. In such 
a conditio n, the gas beco mes locally unstable and develops inhomo¬ 
geneities (IShavivll2001ah that have the effect of reducing the effec¬ 
tive opacity with respect to its microscopic value even before than 
the Ed dington limit is reached. Followine iDotan. Rossi & ShavivI 
(l201lh . we model the effective opacity KcB as: 


— 



r > 0 . 8 , 
r < 0 . 8 , 


(31) 


where F = L,/Lbbb is the local Eddington ratio calculated using 
equation[TT]with M{r), while k is the microscopic opacity: 


k{T) 



1 + (T/To)-i3’ 


(32) 


where Tq = 8000 K. This op acity models the results for pristine 
gas bv iMaver & Duschl dTOOSll . The effective opacity corresponds 
to a an effective Eddington ratio Feff = 1 — 0.16/r when F > 0.8, 
i.e. the gas is effectively sub-Eddington, though it would be super- 
Eddington with the microscopic opacity. Throughout this radiative 
layer (that usually encompasses a tiny fraction of the total mass), 
we assume once again hydrostatic equilibrium (being effectively 
sub-Eddington) and we solve equations [27] and 1281 but we evolve 
the temperature by mean of the radiative gradient with ACefi: 

dP SkcbpL, 

dr IGivacr^T^ 

The luminosity L, remains constant since no energy sources/sinks 
are present within the convective envelope or the radiative layer. 
The inhomogeneities in the radiative layer can maintain the lumi¬ 
nosity sub-Eddington as long as they remain optically thick. Since 
those inhomogeneity have a size of order of the local density scale- 
height, we can estimate their optical depth as t^b ~ XP^esh, 
where x is the ionisation fraction calculated from the Saha equa¬ 
tion assuming equilibrium, and h — |p/(dp/dr)| is the density 
scale-height. Then, the radiative layer extends up to rrad such that 
reff(rrad) = 1- We note that ( typically decreases quickly through¬ 
out the radiative layer, reaching values (^ <C 0.01 at rrad. 

We can finally solve and connect the wind model described in 
Section 12 In particular, we use M* = M(rrad) and i?* = rrad. 
From these quantities we can evaluate {3 associated to the star. As 
in Section 12 we do not model the initial acceleration of mass ex¬ 
plicitly. Instead, we assume that this occurs very quickly around 
rrad, which represent the interface between the hydrostatic part and 
the wind. Then, we assign E = L, lo guarantee the conservation 
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of energy at the interface, because below r^ad there is no net dis¬ 
placement of mass and the total luminosity transported is just L,. 
Finally, we need to specify Mwind- As described in Section [23l we 
assume that the wind connect with the hydrostatic part with a fixed 
Ad* < 1. At the same time, we assume continuity for the density 
and pressure at Trad, which implies: 

Mwind = 4 TT rradP*A4*c*, (34) 

where p* = p(rrad) and c* = Cs(rrad) are the density and the 
sound speed evaluated at rrad as given by the integration of the 
radiative layer, respectively. Once we have /3, E and M, we can in¬ 
tegrate the wind equations as described in Section|2j^and|2j3 with 
the only difference that we use the temperature-dependent opacity 
of equation We evaluate it using as a proxy for the local tem¬ 
perature T = {U which is correct only in the optically-thick 
part of the atmosphere. We discuss the limitations of these assump¬ 
tions in Section[5] 


3.2 Numerical integration 


We proceed to describe the numerical strategy to solve the 
equations in Sec tion 13. 1[ similar to what has been done by 
iDotan. Rossi & ShavivI 1 201 ih . 


(i) We choose the black hole mass M,. 

(ii) We choose one of the central quantity, in particular the total 
central pressure Pc. 

(hi) We need a second quantity to specify all the boundary 
conditions at the centre. Therefore, we guess the value of = 

Pgas.c tP^- 

(iv) We calculate the central quantities: E = 3(1 — (^c)Pc/a, 
Pc = PcCcUmp/iksTc), Ca,c = ((4/3)(1 - Cc) + (5/3)Cc)Pc/pc 
and rn = GM,/ (2Cs^c); we evaluate L, according to eauationl2^ 

(v) We integrate the convective envelope, namely equations 
[^and[^ from the centre (Srs) outward, until we reach Vconv We 
actually integrate the equation in their Lagrangian form, using the 
enclosed mass M (r) as the independent variable. 

(vi) We integrate the radiative layer equations and 1331 

assuming continuity with the convective envelope from rconv till 
frad. We integrate the equations using P as the independent vari¬ 
able since it varies more than the other quantities throughout the 
radiative layer. 

(vii) We calculate the necessary quantities to specify the prop¬ 
erties of the wind using the values at r^ad as discussed in Section 
13.11 namely ft, E and Mwindl specifically, we assume A4* = 0.8. 

(viii) We calculate the wind solution and we check that/whether 
it reaches r^ad self-consistently as described in Section [2(^ 

(ix) We compare the density (rrad) obtained at the end of 
the integration of the radiative layer with the density p*'"‘“‘^^(rrad) 
obtained by the wind integration and we modify in order to 
match the two values. 

We find empirically that the ratio p^"‘"‘^^(rrad)/p*'^'“^^(rrad) 
crosses the value 1 extremely steeply while varying fc and is not 
monotone far from the solution. This occurrence makes difficult 
to use classic methods such as bisection unless the initial guesses 
for the values of /c that bracket the final solution are very close 
to the latter. To overcome this problem, we proceed as follow: we 
choose an initial guess for /c by solving the following equation 
that comes from the scaling r e lations of th e envelope provided by 
iBe gelman. Rossi & Armitara ll2008h and IPotan. Rossi & ShavivI 

( l201lh : 


1.134 (1 - 
mi®/® Cf/® ’ 


(35) 


where Pc = Pc ,7 x 10^ erg cm“® and M, — m, Mq. Then, 
we build a grid of models for several values of /c around the ini¬ 
tial guess, and we progressively refine this grid around the solu¬ 
tion. When we bracket the true solution with a relative precision 
~ 10“®, w e use this br a cketing as the sta rting points for a Brent 
root-finder llBrenllflQT^ : IPress et al.ll2002h . The typical solutions 
of are < 0 . 01 ; indeed, quasi-stars are radiation-pressure domi¬ 
nated in their interiors. Moreover, we note that / typically decreases 
throughout the radiative layer, reaching values f <C 0.01 at Trad. 
This behaviour justifies our simplifying assumption of neglecting 
the gas pressure in the wind, since the ratio 3fcBp/(pmpar®) re¬ 
mains effectively -C 1 through the wind and within the photo¬ 
sphere. Calculations a posteriori of C, in the wind show a decreasing 
behaviour. This suggests that our treatment is at least consistent. 
Of course, the radiation dominated assumption limits our results to 
very massive stars and cannot be extended to e.g. Population III 
stars. 


3.3 Results for quasi-stars 

We run a grid of models exploring a wide range of M, and 
Pc, which maps into M*. Figure [5] summarises our findings in 
the M. — M* plane. Specifically, the various panel shows iso¬ 
contours of the outflow rate Mwind, the accretion onto the black 
hole Mbh = L,/{ric^), where the radiative efficiency 77 ~ 0 . 1 , 
the photospheric luminosity Lphot, and the effective, photospheric 
temperature Tphot. 

Such a plane is characterised by three regions. For high black 
hole masses and relatively lo w envelope masses, no hydrostatic so - 
lution can be found (see also iBeeelman. Rossi & Armita ^l2008h. 
That i s because the quasi-stars would stay beyond the Havashil 
(Il96lh track, which represent a lower limit to the effective tem¬ 
perature (around 4000 K) of a convective envelope in hydrostatic 
equilibrium. Beyond such a limit, no solutions for the hydrostatic 
en velope exist. The second region is the evaporation strip identified 
bv IPotan. Rossi & ShavivI ll 201 ih . This region lays where the evap¬ 
oration timescale fwind = M*/Mwind, i-S- the typical timescale 
over which the stellar envelope would be blown away by the winds, 
is shorter than the accretion timescale t, = M,/Mbh. Within this 
region, the envelope mass is dispersed by the wind before than the 
black holes can accrete further. The third region is the growth re¬ 
gion, where the black hole can accrete substantial mass from the 
envelope (t, < fwind)- 

We find that the iso-contours of Mwin d due to our wind model 
of Sec tion |2] are different from those in IPotan. Rossi & ShavivI 
( I 2 OI ih . Most notably, we do not see any effect of photon-tiring, 
consistently with our result of Section [23] that the ultimate source 
of the wind kinetic luminosity is the advection luminosity due to in¬ 
ternal energy. Specifically, our contours increase mildly when M, 
decreases, while their would become suddenly much more steep at 
the onset of photon-tiring, with M„ind almost independent of M*. 
At constant M,, the different shapes of the iso-contours are such 
that Mwind grows faster at high values of M* until it reaches the 
adiabatic limit discussed in Section l33] 

When our models hit the adiabatic limit, our quasi-star struc¬ 
tures are all characterised by s^, = ,^/iiesc very close to the value 
associated to A4* ~ 0.8 for an adiabatic wind, as shown by Figure 
lAll in the Appendix I aI Above this limit, our choice of a specific 
surface Mach number (« 0.8) over-constrains our mathematical 
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Figure 5. Properties of our quasi-star models in the M, — Af* plane. From upper-left panel, clockwise: thin, blue continuous lines show the isocontours of 
the outflow rate of tho accretion rate on to the black hole M,, of the photospheric bolometric luminosity Lphot. aod of the photospheric effective 

temperature Tpjjot > respectively. The blue, dashed lines in the upper-left panel show the extension in the adiabatic wind regime of some isocontours as discuss 
in the text. The adiabatic wind region is labelled explicitly in this panel, while the label “Growth Region” is omitted for clarity. In each panel, the white, light 
grey and dark grey shaded regions represent the gr owth region, the evaporation strip (i.e. where < f.), and the region where no hydrostatic solutions 

exist as reported bv IPotan. Rossi & ShavivI ilOl ih . respectively; on the other hand, the thick, red, dashed line shows the limit of the evaporation strip (i.e. 
where fwind = f •) for our models. 


system and solutions cannot be found. Technically, we would need 
to leave Af * free to vary and the likely result would be solutions 
with nearly adiabatic winds. 

To test this expectation, we then change method and consider 
the idealised case of a purely adiabatic wind solution (see Ap¬ 
pendix 0 and match it to the hydrostatic envelope. These winds 
are completely specified by the values of the density p(rrad) and 
of the sound speed Ca(rrad) at the end of the envelope. We then 
check that s associated to Ca(rrad) is in the interval 1/6 < s < 
1/4 in order to have a solution with yVl* < 1 (Figure lAn . Fi¬ 
nally, we calculate the adiabatic mass loss Mad according to equa¬ 
tion lAT] and we choose the equilibrium solution (i.e. the value of 


^c) that allows us to match the luminosity L, to the luminosity 
carried by the adiabatic wind evaluated at 7?* = rrad, namely 
3MadCs_*(l + Alad(cs,*) “ ^^ 30/(603 *)). We find a fairly smooth 
transition between the two kinds of models, confirming the nearly 
adiabatic nature of the expected wind when this limit is exceeded 
(see dashed lines in Figure|5] top left panel). 

This difference has also the effect of displacing slightly the in¬ 
terface between the evap oration strip and the growth reg ion, i.e. the 
“threshold-growth line” jPotan, Rossi & ShavivlIroTHl . Our mod¬ 
els predict that such a line moves by a factor of ~ 2 — 3 toward 
higher M,, decreasing the thickness of the evaporation strip. More¬ 
over, for every point (M,,M*) around and within the evapora- 
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Figure 6. Evolutionai'y tracks for quasi-stars. Upper panel: evolution¬ 
ary tracks in the M, — M* plane (see Figure for a description, but 
note that here we plot the evaporation strip determined from our models). 
The green line corresponds to a quasi-stars with initial M, = 100 Mq, 
M* = 10^ 1^0. ^in = 0.1 M 0 yr“^. The blue and red lines both 
correspond to a quasi-star with initial M, = 100 Mq, M* = 10® Mq, 
but with Min = 0.1 Mq yr“^ and Min = 100 Mq yr“^, respectively. 
Lower panel: evolution of M, and M* as a function of time. Continuous 
and dashed lines (coupled with colours) show the evolution of M* and M,, 
respectively. 


ti on strip, we find an Mw i nd sm aller by a factor < 10 compared 
to lPotan, Rossi & ShavivI | |201 ih . which implies a less sudden re¬ 
moval of the envelope when the quasi-stars enters the evaporation 
strip. On the other hand, the iso-contours of Mbh are almost inde¬ 
pendent of M, and they are similar to previous findings. 

The absence of a photon-tired wind has a strong impact on the 
photospheric luminosity of the quasi-stars. Since the wind is mostly 
accelerated at the expense of the internal energy, the diffusive lumi¬ 
nosity coming out at the photosphere Tphot is a large fraction of the 


luminosity L, originally produced by central accretion and trans¬ 
ported through the hydrostatic envelope. Such a luminosity corre¬ 
sponds to Eddington ratios calculated with respect of the whole en¬ 
velope mass that ran ge from Ephot 1 to Ephot 5. 10, even in the 
growth region, where lPotan, Rossi & ShavivI OOl ih find a decrease 
of Fphot due to photon-tiring. Our models predict photospheric lu¬ 
minosities in the interval lO'^® < Tphot/(erg s“^) < lO'^^, with 
iso-contours similar in shape to those of Mwind in the M. — M* 
plane. Such luminosities ar e comparable to moderate bolometric 
luminosities of quasa rs (e.g. lHopkins, Richards & HemauisjlTOOTl : 
iMortlock etZIlioill) and might be observable at high redshift as 
discussed in Section|4]below. At the same time, all our models fall 
in a narrow range of photospheric temperature between ~ 7500 K 
and ~ 9000 K, with temperatures that decrease approaching the 
no-hydrostatic-solution region. We recall from Section l2^ that we 
define the photosphere as the place where the effective tempera¬ 
ture Teff = (i/(47rr^cr))^'^'* equals the proxy for the temperature 
T = {U This happens self-consistently at relatively large 
optical depth r > 10, where T is physically motivated and cor¬ 
rectly influences the optical depth through the opacity law. More¬ 
over, we explicitly check that TeS computed where r = 1 changes 
by at most« 300 K for all our models, suggesting that our determi¬ 
nation of the effective, photospheric temperature Tphot is anyway 
robust. The narrow range of effective temperature is mostly set by 
the microphysical properties of the gas, specifically by the adopted 
opacity law. Indeed, the steep temperature dependence of equation 
I32l is such that the the wind becomes optically thin near the opacity 
drop around Tq = 8000 K. However, the use of the temperature- 
dependent opacity law (equation 1321 in the wind of the quasi-star 
models does not change significantly the main physical properties 
of the wind described in Section |2(i.e. the behaviour of the wind 
when approaching the adiabatic limit), especially within the photo¬ 
sphere where Kes{T) ~ Kes- 

The iso-contours of M„ind, Mbh (or L,), and Lphot behave 
smoothly enough in the M,-M* that they can be reasonably well 
fitted with power-laws. By means of a least-square fitting procedure 
in log-space, we find the following fitting formulas: 

= (1.4 ± 0.1) X 10“'‘ m2 ®® Mq yr“\ (36) 

Mbh = (8.3 ± 0.1) x 10"“ mi "*® m2 °® Mq yr"\ (37) 

Tphot = (3.7 ± 0.1) X 10®® mS ®'^ m2 ®® erg s“\ (38) 


These formulas represent the interpolation between the mod¬ 
els of our grid. Assuming t hat quasi-stars evolve through a se¬ 
quence of equilibrium states i Begelman. Rossi & ArmitaralTOOSl : 


iBall et a ni201ll:lDotan. Rossi & Shavivll201 ih . we can use them to 

calculate evolutionary tracks in the M, — M* plane. We show a 
few example solutions in Figurewhere we assume equations 1361 
and[^to be valid also in the adiabatic regime (only relevant for the 
most massive quasi-star in these examples). Specifically, we solve 
the equations M* = Min — Mvind — Mbh, and M, = Mbh, 
where we allow for smooth constant accretion on to the quasi-star 
envelope through the constant Min. When the quasi-star enters the 
evaporation strip, winds start to become dominant and the envelope 
mass drops while the black hole cannot grow very efficiently. How¬ 
ever, if accretion is intense enough, an almost steady state can be 
established within the evaporation strip, where the accretion from 
outside balances the mass loss due to winds and the black hole can 
grow until the quasi-stars enters the no-hydrostatic-solution region, 
which corresponds to a ratio M*/M, < 20. 
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Figure 7. Predictions of the observability of quasi-stars. Left panel: blue, continuous line and red, dashed line show the probability that a photon with 
wavelength A is captured by the WFC3 camera-nFlbOW filter mounted on HST and by the NIRCam camera+F444W filter planned for JWST, respectively. 
Central panel: predicted AB magnitude in the band of the HST filter F160W as function of the source redshift ^source- Blue-continuous, red-dashed and 
green-dotted lines refer to the photospheric luminosity Lphot = 10^^, 10^®, 10"*® erg s“^, respectively; thin and thick lines refer to the effective temperature 
T’phot = 4000, 8000 K, respectively. Right panel: the same as the central panel, but for the NIRCam camera-HF444W filter planned for JWST. 


4 DETECTABILITY 


We have shown that our models of quasi-stars shine with a bolo- 
metric photospheric luminosity that is bracketed between ~ 10 “*^ 
and ~ 10"*® erg s“^ for M* between ~ 10® and ~ 10^ Mq. At 
the same time, we have found that the interval of photospheric tem¬ 
perature is quite narrow around 8500 K. This information allows 
us to put simple constraints on the detectability of those sources 
by current and future space-based telescopes such as Hubble Space 
Telescope (HST) and James Webb Space Telescope (JWST), re¬ 
spectively. 

Indeed, we can estimate the flux in the filter band X as the 
magnitude: 


mx = Mx + rf( ^source ) + Kx{ •S^source) 5 


(39) 


where Mx is the absolute magnitude in band X, d(z) is the 
distance modulus at the redshift ^ source of the so urce, and 
Kx{zsouTce) is the Jf-correction (e.g. lHogg et alj[2002l) . The dis¬ 
tance modulus is defined as: 


d{z) = 25 + 5 logic ’ (40) 

and it encapsulates all the dependences on the cosmology through 
the luminosity distance Dl(z): 


Dl{z) = il + z)cHf 


Jo 


dz' 


'\/fJm(l + z')^ -f Da 


(41) 


Here and in the following we assume the present day values Ho = 
67.7 km s“^ Mpc~^, flm = 0.309, and Ha = 0.691 for the Hub¬ 
ble parameter, the matter density, and the density of the cosmo¬ 
logical constant, respect ively. These values are consisten t with the 
latest Plank cosmology llPlanck Collaboration et al.ll201.5h . The ab¬ 
solute magnitude is the flux as if the source were 10 pc away from 
the observer: 


Mx = -2.5 log 1 


f 


Lv Txjy) du 
47 r(10 pc)^ V 

V 


p-\-oo 

/ g^Tx{v) 

L Jo 

where is the intrinsic spectral luminosity density (i.e. = 


(42) 


dL/dv) of the source, Tx{v) is the probability of a photon to get 
counted at frequency u with the filter X, and = 3631 Jy is 
the constant spectral flux density of a hypothet ical reference sourc e 
used to ex press magni t udes i n the AB system dOke & Gunnlll983h . 
Following lHogg et al.l(l2002h . we write the JF-correction as: 

Kx{z) = -2.51ogio(l + 2 ) 



The crucial ingredient is the spectra luminosity density of 
the source. We assume that can be modelled as a black body at 
the temperature Tphot, emitting the total luminosity Lphot- Explic¬ 
itly, we have: 

j _ 151/phot {v/vth)^ 

TT^r-th exp(!y/r'th) - 1’ 

where vth = feeTphot/fi. and h is the Plank constant. We discuss 
the limitations of such an assumption in Section|3 

Having an effective temperature ~ 8000 K, our quasi¬ 
star models are expected to be fairly blue; on the other hand, 
cosmologically-motivated calculations pred ict that quasi-stars pop¬ 
ulate mostly massive halos at 2 > 10 dVolonteri & BegelmanI 
l 2010 h . with the consequence of displacing the bulk of quasi-star 
emission in the near-infrared wavelengths > 2 gm. Therefore, we 
focus our analysis on the wide filters in the near-infrared at the 
longest wavelength and contemporary with the highest (effective 
or predicted) sensibility available for HST and JWST. Specifically, 
we consider the filters F160W of the WFC3 camer^f] mounted on 
HST and the filter F444W of the NIRCanfl designed for JWST. 

The left panel of Figure|7]shows the throughput of the consid¬ 
ered filters. The HST data consider also the coupling between the 
camera and the filter, while the JWST data are the predicted trans¬ 
mittance of the filter only. In order to mimic the effect of the cou- 

® http://svo2.cab.inta-csic.es/svo/theory/fps3/ 

® http://www.stsci.edu/jwst/instruments/niream/ 
instrumentdesign/filters/ 
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pling with the camera for JWST as well, we conservatively multiply 
the filter transmission by the fudge factor 0.6, obtaining a maximum 
response similar to the HST values ~ 0.5. 

The central and the left panel of Figure [7] show the predicted 
flux observed by HST and JWST, respectively. We explore the ef¬ 
fect of changing the total luminosity of th e quasi-star and its ef¬ 
fective temperature Tphot to compare with IVolonteri & BegelmarJ 
( l2010h . They assumed an effective temperature of 4000 K, while 
our fiducial model has 8500 K. We compare two different exposure 
times for HST and JWST, namely fexp = 10 h and texp = 10^ s, 
respectively. For comparison, the longest exposure in the F160W 
band of the Hubble Ultra Deep Field ’09 captured with the WFC3 
camera is « 41 hourfl This choice provides a very similar magni¬ 
tude limit of « 27.8 for both instruments and filters. 

We find that hotter quasi-stars are brighter at higher redshift in 
the considered bands. This is because at Tphot = 4000 K, the peak 
of the spectrum is red enough that at high redshift it gets displaced 
beyond the band limit of the filters. This happens at Zsource > 3 
for HST, while at Zsource > 5-6 for JWST, because the F444W fil¬ 
ter extends more in the near-infrared than the F160W one. Quasi¬ 
stars with Tphot > lO'^® erg s“^ could be in principle detected 
by both HST and JWST with the considered integration times at 
^source ^ 10. However, they are close to the magnitude limit for 
HST, while they are well above the same limit for JWST. This sug¬ 
gests that it is fairly unlikely that HST has already observed such 
a source even within a Ultra-Deep-Field-like exposure, whereas 
we expect JWST to be able to detect quasi-stars in the luminos- 
ity range 10"^"^ — 10"^® erg s~ ^ at ^source > 10. Nonetheless, as 
IVolonteri & BegelmanI ll2010li have shown, there might be some 
rare events of quasi-stars forming at redshift as low as Zsource ~ 4. 
Though the bulk of the population is expected to be in place at 
higher redshift, HST might still observe such an outlier. However, 
we caution that those numbers represent the most optimistic esti¬ 
mates since we are neglecting the effect of the environment where 
quasi-stars are expected to live. Indeed, when quasi-stars are har¬ 
boured within gas-rich environment, part of their radiation might 
be absorbed and reprocessed in different wavelengths. 


5 DISCUSSION AND CONCLUSIONS 

Our work addresses the formation of supermassive black hole 
seeds, a major open issue in galaxy formation and high energy 
astrophysics. In particular, we contribute to the assessment of the 
massive seed scenario from direct collapse of gas at the centre of 
(proto)galaxies via the quasi-star mechanism. At this stage in time, 
the assumption that supermassive black holes are grown out of mas¬ 
sive seeds formed at redshift > 10 is theoretically and, of course, 
observationally far from being proved. Whether massive seeds are 
possible, what is their mass function at birth and the prospect of 
detectability with future instruments should be critically assessed 
by a careful investigation of the physical processes at work. 

In this context, we focus on the possible super-Eddington 
growth of embryo black holes inside massive quasi-hydrostatic en¬ 
velopes (quasi-stars), and we constrain the impact of outflows onto 
the final mass with which the black hole seed would emerge at 
the end of this rapid growth p hase. This was already addressed in 
IPotan. Rossi & ShavivI i201 ih but here we consistently solve the 


https : / /archive. stsci . edu/prepds/hudf 09/ 


full equations for a radiation dominated wind, including both dif¬ 
fusion and advection of energy in the flow. We find that such winds 
are ultimately powered by advection luminosity within the flow up 
to the limit where they become nearly adiabatic (i.e. diffusion lu¬ 
minosity is negligible). When applied to quasi-stars, we find that 
the final black hole masses are larger but only by a factor of a few 
with respect to lDotan. Rossi & ShavivI ( 1201 ih : in contrast, the ob¬ 
servable appearance of quasi-stars is expected to be different. They 
are luminous (10“^"^ — 10^^ erg s“^) blue (effective temperature 
of ~ 8000 K) objects. In colour, they differ from predictions that 
ignore ma ss losses, where temperatures can b e up to a factor of 
two lower dSegelman. Rossi & Armita^l2008h . Their characteris¬ 
tics make them promising targets for JWST, while HST may only 
have detected rare, relative closer (z < 10) objects. 

Although our wind treatment improves over previous models 
of quasi-star winds, some caveats need to be discussed. Our steady- 
state assumption implies that we are not modelling the acceleration 
region between the hydrostatic envelop and the wind. It is neces¬ 
sary for us to introduce a discontinuity in the physical quantities, 
and our approach is to assume a vanishingly small acceleration 
scale, where the velocity jumps from u ~ 0 to w ~ 0.8cs. In a 
porous radiation dominated atmosphere, acceleration starts at a ra¬ 
dius where inhomogeneities become optically thin on a scale of 
their size. Since their size is comparable to the density scale-hight, 
it is also reasonable to assume that significant acceleration happens 
on that length-scale. In our simulations, we find that the density 
scale-hight at the onset of the wind is typically 1% of the radius, 
which supports our simplification of an “impulsive” discontinuity 
in velocity. Of course, only time-dependent simulations can defini¬ 
tively prove the correctness of this assumption. We also tested the 
dependence of our models on our working value of A4* = 0.8. 
As long as A4* > 0.5, the main results are weakly sensitive to the 
exact choice (see also Figure and Section |2. 2. It . However, for 
lower values of AI*, solutions can still be found, but show unphys¬ 
ical behaviours when approaching the adiabatic limit. Consistently 
with the preference for A4^, > 0.5 shown in Section 12.2.11 this 
suggests that a value Al* 1 necessary to better satisfy the 
assumption of a steady state from the base of the wind. 

In contrast with velocity, temperature and density are contin¬ 
uous across the two regions (hydrostatic envelope and wind). This 
allows us to calculate M„ind through the continuity equation. We 
do not adopt the mass loss rate assumed in lPotan, Rossi & ShavivI 
1 20 njubecause it seems to require the presence of a critical point 
I Shavivl2001al f3l. that is instead absent when solving the equations 
of a purely radiation dominated wind (eguations IlISl l. Finally, fur¬ 
ther progress should include the effect of rotation, which in this ver¬ 
sion is omitted, and may instead cause larger mass loss rate around 
the polar axis. However, a funnel around a polar axis may be carved 
by the presence of a jet, that would result in a smaller amount of 
energy being injected into the quasi-star envelope. This last would 
simply have the effect of changing the black hole mass to envelope 
mass ratio of a given solution, since for a given quasi-star mass, 
more accretion energy is needed to support it. More sophisticated 
3D simulations, possibly coupling hydrodynamics with radiative 
transfer, are necessary to validate the overall picture and to assess 
the impact and interplay of such processes; this opens many possi¬ 
bilities for future additional investigations. 

Our assumption of an evolution through equilibria breaks 
down when the quasi-stars encounter the no-hydrostatic-solution 
region. Since our approach does not allow us to infer quantita¬ 
tive predictions of the final evolutionary phase, w e can only specu - 
late about what might happen afterward (see also lSall et alj|201lh . 
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When a quasi-star reaches this region, the envelope is still a fac¬ 
tor of ~ 20 larger than the black hole. Such an envelope can¬ 
not remain in hydrostatic equilibrium and it is conceivable that at 
least a fraction of it may collapse onto the central black hole, lead¬ 
ing to further accretion. Since the infalling gas is expected to be 
radiation-dominated and optically-thick, the associated accretion 
episode might involve a large f raction of mass, at least where the 
black hole potential dominates iBeeelmanll 19781) . This line of rea¬ 
soning suggests that the masses that we infer in Section [53] might 
be lower limits. However, if the infall proceeds out of equilibrium 
on a few dynamical timescales, the inflow rates may be much larger 
than the Eddington li mit of the black hole and feedback may limit 
further accretion (e.e. [Johnson et al]|201lh . 

Finally, to estimate the detectability of quasi-stars, we assume 
a black body spectrum, that may be accurate for our broad band lu¬ 
minosity estimates of this pristine object. Of course, reliable spec¬ 
tral predictions must instead account for lines and electron scatter - 
ing, but this should required a dedicated study te.g. lSchaereJ2002h . 
We also warn the reader that our predictions are for the intrinsic lu¬ 
minosity of quasi-stars and that absorption is not accounted for. 

Despite these caveats, our work confirms that super-Eddington 
accretion onto newly bom black holes within quasi-stars is likely 
responsible for vigorous mass loss, which in turn may limit the 
growth of the black holes. If this is the case, forming massive 
(> 10 "^ Mq) seeds via the quasi-star scenario is therefore more 
difficult than is generally wished for. It requires massive inflows 
of gas at the centre 10 Mq yr“^, see e.g. our Figure]^, that 
are generally associated with massive (> 10® Mq) and therefore 
relative rare halos at high redshift. This is at least qualitatively in 
agreement with recent observational studies that fail in finding can¬ 
didates of massive quasars at 2 > 5 and possibly put constr aints 
on the early black hole formation modes l lWeigel et alj|201^ . To 
assess how rare those massive seeds are, and whether they can still 
account for the bright quasars observed at 2 > 6, a consistent cos¬ 
mological evolution needs to be computed. This is the topic of a 
follow up paper, that will also allow us to make a more quantitative 
assessment of the detectability of these fascinating objects. 
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APPENDIX A: REVIEW OF ADIABATIC WINDS 


For convenience, we briefly review the main features of an adi¬ 
abatic wind characterised by the equation of state p = , 

whe re AT is a constant (for further and more general readings, see 
e.g. lUolzer & Axfordll97(]li . Within the formalism introduced in 
Section im the equations of such an isentropic wind are: 


l--)w' = -l-b 
w / 



(Al) 


w' &s 1 = 0. 


(A2) 


These equations can be derived by our starting equations of Section 
I2.1l in the optically-thick limit when k —>■ +oo and L{r) —>■ 0. 
From the equation above we immediately see that the critical (or 
sonic) point (i.e. where w = s), when present, coincides with the 
singular point where w' can diverge. To avoid that, we require that 
the left-hand side of eouation lAll is 0 at the critical point Xc, which 
implies the relation between Xc and the critical velocity Wc'. 

Xc = l— 4uic. (A3) 


This equation shows that the maximum speed to have a critical 
point outside the stellar surface (i.e. Xc > 0, which corresponds to 
a subsonic solution at i?*) is Wc 1/4, or ^ GM*/(2i?*). On 
the other hand, a critical point approaching infinity (i.e. Xc 1) 
corresponds to a critical velocity Wc —>■ 0. 

EQuation lA2l describes the conservation of energy and entropy. 
The request of a non-diverging critical point sets the total energy 
associated to wind: 


e = 'u;-l-6s-l-x = StXc -I- 1. 


(A4) 



Figure Al. Relation Al* ad('5*) between the sound speed and the Mach 
number at 7J* for an adiabatic wind (thick continuous line). The thick, 
dashed line shows the scaling of the mass outflow rate M^d with the sound 
speed Cs^*. The thin, dotted line shows the position s* = 1/6, where the 
Mach number approaches asymptotically zero. The thin, dot-dashed lines 
show that when s* = 1/4, the wind has AJ* = 1 exactly at R*. 


This condition, when evaluated at R ^,, provides a relation between 
Wc and the quantities w* and s*. The second condition that al¬ 
lows a full determination of w* and s* given Wc is the conserva¬ 
tion of mass, which reads as equation [T] combined with the rela¬ 
tion p/p* = (s/s*)^, which comes from the equation of state and 
the general definition of sound speed in equation [25] The final and 
complete relations between Wc and ui* and s* are: 


ui* -b 6s* = 3wc -b 1, (A5) 

16 wy^s* = (A6) 

When Wc —> 0, the second relation shows that w* s* —>■ 0 and 
either w* or s* has to go to zero. Since this limit case corresponds 
to have a critical point well outside Jf*, we have w* <C s* and 
therefore ui* —>■ 0. The first relation then imply that every adia¬ 
batic solution with a critical point at a finite radius has necessarily 
s* > 1/6, while every solution with a critical point at a finite radius 
larger than i?* requires 1/6 < s* < 1/4. Combining equations 
IA5l and lA6l provides also a unique relation AJ*,ad(s*) between the 
sound speed and the Mach number at 77*. This relation is shown in 
Figure IaTI 

We can finally estimate the mass outflow rate through the 
conditions at the critical point and map them back to the prop¬ 
erties of the flow at Jf*. At the critical point, the outflow rate 
reads M = 47rrcp(rc)uc. Using the relations Xc = GMi,/{2Vc), 
p(xc) = p*(i’c/cs,*)® and equation |A5 1 in dimensional form, we 
obtain: 


Ma.d = 


47rG^M*V* 

ysc®* 


V 6 



(A7) 


This depends on the properties of the star, such as M* and Vesc, on 
the density p* that sets the normalisation, and on the sound speed 
Cs,*. Figure lATI shows also the scaling of Mad with the sound speed 

Cs,*. 
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